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Energy spectra, electron densities, pair correlation func- 
tions and heat capacity of a quantum-dot lithium in zero 
external magnetic field (a system of three interacting two- 
dimensional electrons in a parabolic confinement potential) 
are studied using the exact diagonalization approach. A par- 
ticular attention is given to a Fermi-liquid - Wigner- solid 
transition in the ground state of the dot, induced by the intra- 
dot Coulomb interaction. 

PACS numbers: 73.21.La, 73.20.Qt, 71.10.-w 



I. INTRODUCTION 

Quantum dots0 are artificial electron systems (ES) re- 
alizable in modern semiconductor structures. In these 
systems two-dimensional (2D) electrons move in the 
plane z — in a lateral confinement potential V(x,y). 
The typical length scale Iq of the lateral confinement 
is usually larger than or comparable with the effec- 
tive Bohr radius of the host semiconductor. The 
relative strength of the electron-electron and electron- 
confinement interaction, given by the ratio A = Zo/ a Bj 
can be varied, even experimentally, in a wide range, so 
that the dots are used to be treated as artificial atoms 
with tunable physical properties. Experimentally, quan- 
tum dots were intensively studied in recent years, usir 
a variety of different techniques, including capacitance 
transports, far-infraredo and Raman spectroscopy!! 

From the theoretical point of view, quantum 
dots are ideal physical objects for studying effects 
of electron-electron correlations. Different thec 
ical approaches, including analytical calculation^ 
exact ^a^nalizatior£3c3, quantum_Monte Carlo 
(QMC)ci^L density functional theoryErE] and other 
methodsE3c3, were applied to study their properties, for 
a recent review see Ref.LJ. Until recently most theoreti- 
cal work was performed in the regime of strong magnetic 
fields, when all electron spins are fully polarized. In the 
past three years a growing interest is observed in study- 
ing the quantum .dot-properties in zero magnetic field 
B = ome-MNHM The aim of these studies is to 
investigate the Fermi liquid - Wigner solid crossover in 
the dots, at a varying strength of Coulomb interaction. 
Detailed knowledge of the physics of such a crossover 
in microscopic dots could be— compared with that ob- 
tained for macroscopic 2DESc3 , CJ and might shed light 
on the nature of the metal-insulator transition in two 



dimensionsEl 

So far, full energy spectra of an ./V-electron parabolic 
quantum dot in zero magnetic field, as a function of the 
interaction parameter—^, were published only for N = 2 
(quantum-dot heliumtH). For larger N a number of re- 
sults for the ground state energy of the dots were reported 
at separate points of the A axis. Sometimes, however, re- 
sults obtained by different methods contradict to each 
other, and full understanding of physical properties of 
iV-electron dots at B = has not yet been achieved. 

Even for three electrons in a parabolic confining po- 
tential (quantum-dot lithium) available in the literature 
results are somewhat confusing and do not give a clear 
and exact picture pL the ground state of the dot. Egger 
et al. reported inEIJ on a transition from a partly to a 
fully spin-polarized ground state at A w 5 caused by the 
intra-dot Coulomb interaction (z-component of the total 
spin 5*°* changes from 1/2 to 3/2). This conclusion was 
based on results of QMC calculations (multilevel block- 
ing algorithm). For the 1/2 state calculations gave at 
A = 4 slightly lower energy, than for the 3/2 state, and 
higher energies at A-j= 6, 8 and 10. In a subsequent (Er- 
ratum) publication^ other (corrected) values for the en- 
ergy of these states were reported. At A = 4 they gave the 
same energies of the 1/2 and 3/2 states. The conclusion 
from these two publications in the part concerning the 
three-electron quantum dot is thus not clear. HauslerQ, 
calculating the energy spectrum of a three-electron dot 
using the so-called pocket state method, gave a more def- 
inite indication on the existence of such a transition. He 
presented his results, however, in some specific for the 
method and not-directly related to A units, which does 
not allow one to quantitatively characterize the transi- 
tion (e.g. to get the value of the interaction parameter 
A at the transition point). Hausler also pointed out that 
some indications on such a transition can be also seen in 
earlier studies of a three-electron dot: in Ref.EJ (A 2) 
the ground state of the dot at = was found to be 
partly spin-polarized, while in Rcf.cil with a substantially 
larger value of A the fully spin-polarized state turned out 
to have the lower energy. As seen from this brief out- 
line, available data indicate that this transition seems to 
exist but it is not yet quantitatively understood, and its 
physical origin is not completely clear. It is also not clear 
whether and how this transition in the spin state of the 
system is related to the Fermi-liquid - Wigner-molecule 
transition in the dot. 

In this paper I present results of a complete theoretical 
study of a three-electron parabolic quantum dot. Using 
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exact diagonalization technique, I calculate full energy 
spectrum of the dot, as a function of the interaction pa- 
rameter in the range < A < 10. At A = A c = 4.343 I 
find a transition in the ground state of the dot, accompa- 
nied by the change of the total spin quantum number. I 
study the densities and the pair-correlation functions in 
the ground and the first excited states of the system at a 
number of A-points, including the vicinity of the transi- 
tion. These results show that physical properties of the 
dot dramatically change at the transition point, corre- 
sponding to properties of Fermi liquid at A < A c , and of 
a Wigner molecule at A > A c . I also calculate some ther- 
modynamic properties of the dot: the heat capacity and 
the volume-pressure diagram. The temperature depen- 
dence of the heat capacity clearly exhibits characteristic 
features, related to the transition, in the temperature 
range corresponding to a few K for typical parameters of 
GaAs dots. Other experimental consequences from pre- 
dictions of this paper are also discussed. 

In Section [n I briefly describe the model and the 
method used in calculations. Results of the work are pre- 
sented in Section H. Concluding remarks can be found 
in Section IV. 



II. MODEL AND METHOD 

A. The Hamiltonian 

I consider three 2D electrons moving in the plane z = 
in a parabolic confining potential V(r) — m*u>Qr 2 /2, 
r = (x,y). The Hamiltonian of the system 
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(N = 3) commutes with operators of the total angular 
momentum L' ot , (squared) total spin Sj ot and projection 
of the total spin S^ ot on some (£-) axis (not necessarily 
coinciding with the z-axis). This gives three conserving 
quantum numbers L to t = £| ot , S to t and S^ ot . No mag- 
netic field is assumed to be applied to the system. 



and the spin functions Xs{&)- Here Z = ^Ti/m*uj is 
the oscillator length, and (n, I, s) are the radial, asimutal 
(angular momentum) and spin quantum numbers of the 
single-particle problem (n > and I are integer, s — 
±1/2). All the single-particle states (|J) can be ordered 
and enumerated, e.g. 0i = (nls)i — (0, 0, |), 02 — (0, 0, I 
), 03 = (0,1, |), etc. The energy of the states @ does 
not depend on spins, 

E nls = htu (2n+ \l\ + 1). (4) 



C. Basis set of many- particle states 

A complete set of many-particle states ^f u , u = 1, 2, . . . 
is formed by placing particles in different single-particle 
states, e.g. *i = |0i0 2 <^3>, *2 = 1010204), • where 
\4>a<Pp<P-/} are Slater determinants 
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and £j = (r^, <jj). All the many-particle states ^f u can be 
also arranged, e.g. in order of increasing of their total 
single-particle energy 



E°p = E, 
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(6) 



and enumerated. 



D. Solution of the Schrodinger equation 

Expanding the many-body wave function in a complete 
set of many-particle states, 



(7) 



I get the Schrodinger equation in the matrix form, 

J2( H uu> -E5 UU ,)C U , =0. (8) 



B. Basis set of single- particle states 

A complete set of single-particle solutions of the prob- 
lem 

0„ is (r, a) = (f nl (r)xs (c) (2) 
is the product of the Fock-Darwin orbitalsEU 

1*1 
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The conservation of the total angular momentum Ltot 
and the projection of the total spin S^ ot allows one to 
chose the many-body states for the expansion (^) under 
additional constraints 
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h — Ltot, Si — S t 



tot 



(9) 



This reduces the size of the matrix in (||) and facilitates 
calculations. 

Numerically diagonalizing the eigenvalue problem (^) 
I get a set of energy levels 
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and corresponding eigenfunctions 



L toU S\ ot ,m 



(10) 



(11) 



as a function of the interaction parameter A = Ig/ag- 
The number m — 1,2,... enumerates the energy levels 
of the system in the subspace of levels with given L to t 
and S(- ot . After the diagonalization problem is solved, 
the eigenvalues of the total spin are calculated for each 
level m from 

S tot (Stot + 1) = (*i tot ,S*",m|St ot |*Z tot ,S*<-,m)- (I 2 ) 

All the matrix elements of the Hamiltonian H uu i and of 
the operator (S 2 ot ) uu / are calculated analytically. 

All the energy levels with non-zero L to t and S*^ ot are 
degenerate, 
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at = E. 
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(13) 

Presenting below results for the energy of the states 
(Ltot, Sj ot ), I omit the corresponding signs [for instance, 
(1, 1/2) stands for (±1, ±1/2) with all possible combina- 
tions of signs]. Degeneracy of levels are calculated ac- 
counting for (O). 

E. Properties of the states and the heat capacity 

After the Schrodinger problem is solved and all the 
energy levels and the eigenfunctions are found, I calculate 
the density of spin-up and spin-down polarized electrons 
in the states (L to t, S^ ot ,m), and the corresponding pair- 
correlation functions. These quantities are calculated as 
averages of the operators 
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and 
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P ffCT ,(r,r') = ^ ^ <5(r-r 4 )<5(r'-r,)<W^. (15) 

with the eigenfunctions ^L tot! s* ot ,m- All the matrix ele- 
ments of the operators (|l4|) and ([l5]) are calculated ana- 
lytically. 

As the method offers an opportunity to find all the 
energy levels of the system, one can also calculate ther- 
modynamic properties of the dots. I calculate the heat 
capacity as C\ — (dE/dT)\, where 

J2 n En(\)e~ E ^/ T 



E = E(X,T) = 



(16) 



T is the temperature, and the sum is taken over all (low- 
lying) energy levels accounting for their degeneracies. 



F. Convergency of the method 

The number of all many-particle states in the problem 
is infinite, and the size of the matrix H uu i in Eq. (||) is 
infinite too. To perform practical calculations I restrict 
the number of many-particle states in the expansion (0) 
so that the total single-particle energy (||) of the involved 
many-body states is smaller than some threshold value 
E t h, E^ p < E t h- The larger the threshold energy E t h, the 
broader the range of A in which results are convergent and 
reliable. Typically, less than 1000 many-particle states 
were sufficient for all the calculations presented below. 

Convergency of the method is illustrated on Figure ||, 
where the energy -£(1,1/2) of the lowest state (to = 1) 
with (Ltot,Stot) = (1, 1/2) is shown as a function of 
A = lo/as for increasing threshold energy E t h- The 
curves are labeled by E t h and the number of many-body 
quantum states N m i, s involved in the expansion (|7|) . One 
sees that including about 1000 many-body states leads 
to very accurate results for the energy at A < 20. Notice 
that below I present results for the energy in the interval 
A < 10, where the method is practically exact: at A = 10 
I found that E {1a/2) /Tlujo = 17.627891 at N mbs = 1024, 
and 17.627974 at iV m f, s = 549. The difference comprises 
4.7-10~ 4 %. 



III. RESULTS 

All the lengths in this Section are measured in units 
Iq, all the energies - in units ?ioJo, the densities and the 
pair-correlation functions - in units (7rZ§) — 1 and (ttIq) 
respectively. 



A. Energy spectra 

1. Exact results 
The interaction parameter in the problem 
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characterizes the relative strength of classical Coulomb 
(~ e 2 ) and quantization effects (~ K). The limit of small 
A corresponds to a weakly interacting system (e 2 — * 0) 
and can be treated exactly. The ground state in this limit 
is realized in the configuration [(0, 0, t)(0, 0, |)(0, 1, f)], 
i.e. (L tot , Stat) = (1, 1/2), with the energy 



lim Ecs/fi^o = 4. 



(18) 



The opposite case A — > 00 corresponds to the purely 
classical limit (U — * 0). In the classical ground state, jelec.- 
trons occupy the corners of an equilateral trianglePlri, 
with the distance 
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Rci = ld/3 1/6 (19) 

from the origin. The ground state energy at A = oo is 

Egs = 3 5/3 e c //2. (20) 

Here l c \ — (e 2 /m*^) 1 / 3 and e ; = e 2 /l c i are the classical 
length and energy units. So, in the classical limit 

lira Eas/hoJo = 3 5/3 A 2/3 /2. (21) 

A — >OG 

Equations (jl|) and ( pl| ) give asymptotes of the ground 
state energy at very small and very large A. 

At arbitrary A one needs numerical calculations. Fig- 
ure |^ shows the low- lying energy levels of a three-electron 
parabolic quantum dot with the total angular momen- 
tum L to t = 1. Shown are only the states which start 
from E/huiQ < 6 at A = 0. The lowest-energy state (the 
ground state in the subspace of levels with Ltot = 1) cor- 
responds (in the limit A — > 0) to the configuration [(0, 0, j 
)(0,0,|)(0, 1,T)], and has the total spin S to t = 1/2. This 
state is four-fold degenerate (Ltot = ±1/2, S^ ot = ±1/2). 

Figure || shows the low-lying energy levels with the 
total angular momentum Ltot = 0. Shown are only 
the states which start from E/Tlloq < 7 at A = 0. 
The lowest-energy state has the total spin Stot = 3/2 
and four-fold degeneracy (S^ = ±3/2, ±1/2). In the 
limit A — * the state with the full spin polarization 
{Stot — 3/2, S^ ot = +3/2) corresponds to the configura- 
tion [(0, -l,t)(0,0,t)(0, 1,T)]- The state with a partial 
spin polarization {S to t = 3/2, S^ ot = +1/2) corresponds 

in the limit A — * to the configuration (*i+* 2 +*3) /\/3, 
where 

*i = [(o,-i,|)(o,o,T)(o,i,T)], 

vF 2 = [(0,-l,T)(0,0,i)(0,l,T)], (22) 

* 3 = [(o,-i,T)(o,o,T)(o,i,i)]. 

I performed similar calculations for L to t from to 9 and 
for all total-spin states Stot = 3/2 and 1/2. Results for 
Ltot from to 2 are shown in Figure |4j. Only the lowest- 
energy levels {m = 1) for each pair of numbers {L to t, Stot) 
are shown in the Figure [this means that, say, a (1,1/2)- 
level with m > 1 (not shown in the Figure) may have 
the lower energy than the exhibited level (2,3/2) with 
rn = 1]. At some critical value of A (A = A c = 4.343) 
one observes a crossing of the two lowest-energy levels 
(1, 1/2) and (0,3/2) (more clearly seen in Figure g where 
the energy difference £(0,3/2) — £(1,1/2) is plotted versus 
A). At the critical point the total spin of the system in 
the ground state changes from S to t = 1/2 at A < A c to 
Stot — 3/2 at A > A c . Near the critical point the gap 
between the ground and the first excited states behaves 
as 

1^(0,3/2) - £(i,i/2)|/^o = 0.02766| A - A c |, (23) 

with a jump of the derivative of the ground state energy 
with respect to A. At large A (more exactly, at A = 10) 
the energy difference is £(0,3/2) —£(1,1/2) ~ —0.04167^0- 



In Table g exact results for the energies of the states 
(1,1/2) and (0,3/2) are compared with QMC results from 
Ref.E3. One sees that the accuracy of the QMC resultsEHl 
is in general very high, but the errors are not always small 
(see e.g. A = 10) compared to the difference between the 
energies of the ground and excited states. 

2. Approximations 

The energy of the states (1,1/2) and (0,3/2) can be 
approximated, in the interval < A < 10, by the formula 

E = E$ s + e cl (AX + \/£ 2 + C 2 A 2 - B) , (24) 

where Eq S is the classical ground state energy ([20]) , X = 
h~u)o/e c i, and the numbers A, B, C for the two considered 
states are 

(A, B. C)(i,i /2 ) = (3.11536, 2.93076, 0.917954), (25) 

{A, B, C) (0,3/2) = (2.84171, 2.44529, 2.13633). (26) 

For the state (1,1/2) the difference between the exact 
energy and the approximation (p4|)-(p5"|) is about 0.83% 
at A = 0, does not exceed 0.22% at 0.05 < A < 10, and 
tends to zero with growing A. For the state (0,3/2) the 
difference between the exact solution and the approxi- 
mation is smaller than 0.44% in the whole range of A. 
It should be noted however, that the difference between 
the energies of the two states £(0,3/2) — £(1,1X2) is re- 
produced by the approximate formulas (|2~i|)-(|26|) with a 
substantially worse accuracy. 

B. Heat capacity 

Figure |^ exhibits the calculated low-temperature heat 
capacity C\ as a function of T and A. About 30 lowest- 
energy levels for each Ltot from to 9, with correspond- 
ing degeneracies, were taken into account in this calcu- 
lation. A pronounced peak related to the gap between 
the ground and the first excited state can be clearly seen 
in this Figure. The value of the peak temperature T p {\) 
depends on A as |£(o,3/2) — £(1,1/2) I (compare to Figure 
H), and disappears at the critical point A = A c . The most 
dramatic variations of the heat capacity are the case in 
the range kT < 0.1?iwo, which corresponds, at a typical 
confinement of GaAs quantum dots {TiUq ~ 3 meV) , to a 
few-K temperature scale. 

C. Electron density and correlations 

Due to rotational symmetry of the Hamiltonian dl]) 

(L S ) 

the density n K a tot ' tot '{r,6) of spin-up and spin-down 
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polarized electrons in the quantum-mechanical states 
(L tot, S tot) does not depend on the angular coordinate 
0, and are shown below as a function of r only. The pair- 
correlation functions P aa ' (r, r') are plotted as a function 
of r/lo, for all orientations of spins, at the second coordi- 
nate r' being fixed at the classical distance (|l^) from the 
origin, r' = (0, — Rd) (the second subscript corresponds 
to the spin of a fixed electron) . 



1. The state (L tot , S to t) = (0, 3/2) 



The total density of electrons n ,(°' 3 / 2 ) (r) = n|°' 3/2) (r) + 
n (°,3/2)^ in the state (L totl Stot) — (0, 3/2) at a few val- 







ues of the interaction parameter A is shown in Figure 
The function n' ' 3 / 2 ) (r) also determines the densi- 
ties of spin-up and spin-down polarized electrons: in the 
state with the total spin projection S^- ot — +3/2 the 

density of spin-up electrons coincides with n(°' 3 / 2 )(r), 
while the density of spin-down electrons is zero; in the 

5* ot = +1/2 state one has n^ 3/2 ' +1/2) (r) = §n(°> 3 / 2 )(r), 

and nf' 3/2 ' +1/2 \r) = ^^^(r). One sees that at all 
A maxima of the electron density lie at a finite distance 
from the origin. At A > 1 they are very close to the 
classical radii (|l9|), shown by triangles on the Figure; at 
A <C 1 (weak Coulomb repulsion) they are at even larger 
r, due to the exchange repulsion. 

Figure |^ exhibits the pair-correlation function 
P TT (r,r') in the state {L toU S to t, Sf 1 ) = (0,3/2, +3/2) 
(three other functions in this state obviously vanish, 
P n = P n = P u = 0; in the state with = +1/2 one 

has P^i\ 1/2) = iP T ( T °' 3/2 ' +3/2) , and P^ /W) = 0). 
The interaction parameter A assumes the values 0.1, 2, 4, 
and 8, from up to down. At small A electron-electron in- 
teraction is weak, and the pair-correlation function has a 
form of a single peak centered opposite to the fixed elec- 
tron. With growing A this peak is splitted onto two ones, 
and this splitting becomes very pronounced at strong 
Coulomb interaction (A > 4). 

The internal structure of the state (0,3/2) thus cor- 
responds to the angle-averaged classical configuration of 
an equilateral triangle (see the inset to Figure ||). This 
structure is highly symmetric and "rigid" : the ratio of 
the sides of the triangle remains unchanged when the cur- 
vature K of the confiment potential varies (K cx A~ 4 ). 
This state is a quantum-dot analog of the Wigner solid 
(a Wigner molecule). 



2. The state (L tot , S to t) = (1, 1/2) 

In the state (1,1/2) the spatial distribution of elec- 
trons is less trivial. Now we have two sorts of particles: 
two electrons are polarized up, and one electron is po- 
larized down. Figure ^ exhibits the total electron den- 
sity (spin-up plus spin-down). One sees that at small A 



electrons behave as non- or weakly interacting particles 
(a Fermi-liquid- type state), forming the structure with 
maximum of the electron density at r = 0. Such a pic- 
ture is the case up to A ~ 2, when a weak minimum 
of the density at r = appears. At even larger A the 
influence of electron-electron interaction becomes more 
important: the density of electrons qualitatively looks 
like in the state (0,3/2), with a minimum at r = and 
maxima close to the classical radii (fl9|). 

Additional and even more interesting information can 
be extracted from Figures 10 and II], which show sepa- 
rately the densities of spin-up and spin-down polarized 
electrons. One sees that at small A the one spin-down 
electron occupies the center of the dot, while the two 
spin-up electrons rotate around the center with a maxi- 
mum of the density at a finite distance from the origin. 
Such a situation is a peculiar quantum-mechanical fea- 
ture: it is not encountered in the classical picture. With 
the growing A (electron-electron interaction increases) 
the two spin-up electrons push the spin-down electron 
out from the center, but the structure "one (spin-down) 
electron is essentially closer to the center, two (spin-up) 
electrons rotate around" conserves up to A ~ 2: the den- 
sity of the spin-up electrons has a clear maximum at a 
finite distance from the origin, while the density of spin- 
down electron still maximal at r = 0. At even larger A 
(> 4) (stronger electron-electron interaction) the density 
maximum of the spin-down electron is shifted to a finite 
r, but at any A it is closer to the origin, than the density 
maximum of the spin- up electrons. 

These features can be also seen in Figures |l2], |l3|, and 
|lj, which exhibit the pair-correlation functions P^, P^, 
and P n in the state (L tot , S tot , Sf*) = (1,1/2, +1/2) 
(-Pu = in this state). Compare for instance Figures [l| 
(Pjl) and [l3] (-Pj.t). In both cases a spin-up electron is 
fixed at the classical distance from the origin. Let A = 2 
(the second row of plots). One sees that the maximum of 
the pair-correlation function is about two times closer to 
the center of the dot for the spin-down electron (Figure 
|l3| ) than for the spin- up electron (Figure ]l^) . At A = 4 
this difference is less pronounced but can also be seen. In 
Figure |l3| one also sees that in the limit of weak Coulomb 
interaction A <C 1 the spin-down electron is localized in 
the center, in agreement with the above analysis of the 
density plots. 

The internal structure of the state (1, 1/2) thus corre- 
sponds to an angle-averaged configuration of an isosceles 
triangle, with two spin-up electrons at the base corners 
and one spin-down electron at the top of the triangle. 
This structure is less symmetric than that of the (0, 3/2) 
state and "soft": the ratio of the sides varies with A, 
changing from 1/2 at A = to 1 at A = 00 (see the in- 
sets to Figure ||). This state is of a Fermi-liquid type 
(Fermi-gas at A — ► 0). 
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3. Reconstruction of the ground state 

Now, we can understand the physical reason of the 
transition (1,1/2) <-* (0,3/2) at the varying parameter 
A. Consider what happens with the ground state of the 
system, when the curvature K ex o>g oc A -4 of the con- 
finement potential varies from small (the limit of strong 
Coulomb interaction A — > oo) to large values (the weak 
Coulomb interaction regime). At small K the quantiza- 
tion effects are negligible, %oj /t c i = A~ 2 / 3 <C 1, electrons 
are at a large (compared to ^o) distance from each other, 
and form a quasi-classical equilateral-triangle structure. 
In this highly-symmetric structure all electrons should 
be equivalent (have the same spin), therefore the total 
spin of the dot in this limit is S to t = 3/2. Increasing 
the curvature reduces the area of the triangle. Its form 
however first remains unchanged. Further increase of the 
curvature costs energy, and the system is reconstructed, 
at A = A c , to another ground state with a more compact, 
"soft" isosceles-triangle structure. In this less symmet- 
ric structure one electron should differ from two others, 
therefore the transition to the new ground state is ac- 
companied by the change of the total spin Stot — * 1/2. 
Further increase of the curvature changes the ratio of 
sides of the isosceles triangle, but not its symmetry. 

Exactly at the transition point A = A c the physical 
properties of the dot change very sharply. Figure |l5| 
shows the total density of electrons n e (r) at A — A c in the 
states (0, 3/2) and (1, 1/2). When the system passes from 
the (0,3/2) to the (1,1/2) state, electrons are pushed to- 
wards the center of the dot with a 50% increase of the 
density at the point r = 0. Figure |l^ shows the area of 
the dot 



A= — 

N 



If 7T 

- / dvirr 2 n e (r) = — (i& G8 \ £j r?|tf G S> 



(27) 



[^GS is th e ground-state wave function] versus the cur- 
vature K oc A~ 4 near the critical point A = A c . As the 
curvature of the confinement potential can be treated as a 
"pressure" acting on electrons of the dot from the confin- 
ing potential, Figure |l6| can be considered as a "volume" - 
"pressure" diagram. One sees that increasing the pres- 
sure leads to a discontinuous jump of the volume (with 
a (5-like peculiarity in the compressibility) at the critical 
point A = A c . In a real system the transition shown in 
Figure 16 may happen with a hysteresis. 



in 2DES; here r s = 1/ \jim s <iB- In such coordinates the 
ground-state-energy curves, both in dots and in 2DES, 
look very similarly: the Fermi-liquid and the Wigner- 
solid energies are smooth, almost linear functions [for 
dots see Eq. (p4l)], intersecting each other at some value 
of the interaction parameter. In the regime of strong 
interaction the fully spin-polarized Wigner solid is the 
ground state in both the 2DES and the dot; in the regime 
of weak interaction the unpolarized (partly-polarized in 
three-electron dots) Fermi-liquid state is the ground state 
of both systems. The transition (0, 3/2) «-> (1, 1/2) in the 
quantum-dot lithium is thus an analog (or a precursor) 
of theJfermi-liquid-Wigner-solid transition in an infinite 
2DESH 

The transition (1, 1/2) <-> (0, 3/2) should be observable 
in many experiments. It should manifest itself in any 
thermodynamic quantity. The difference in the ground- 
state total spin Stot should be also seen in the Zeeman 
and spin splitting of levels in magnetic field, both parallel 
and perpendicular to the plane of 2DES, as well as in 
Kondo-effect experiments. The structure of levels could 
be also studied by Raman spectroscopy. 

The method of the paper can be used for studying sys- 
tems with more particles and/or in a non-parabolic con- 
finement potential^. It can also be used for investigating 
other properties of the system, for instance, the influence 
of impurities, or response of the dot to external fields. It 
is seen already now, for instance, that an asymmetrically 
located impurity will qualitatively differently affect the 
ground state of the dot at A < A c and at A > A c : in the 
former (the latter) case the ground state is degenerate 
(non-degenerate) with respect to L to t, and the impurity 
will result in a splitting (a shift) of the ground state level. 

To summarize, I have performed a complete theoretical 
study of a system of three Coulomb-interacting electrons 
in a parabolic confining potential, and investigated in 
detail physical properties and the origin of the Fermi liq- 
uid - Wigner solid transition in the ground state of the 
three-electron parabolic quantum dot. 



ACKNOWLEDGMENTS 

The work was supported by the Sonderforschungsbere- 
ich SFB 484, University of Augsburg. I thank Klaus 
Ziegler, Ulrich Eckern, Vladimir Sablikov, Teun Klap- 
wijk and Miodrag Kulic for useful discussions, as well as 
referees of the paper for useful comments. 



IV. CONCLUDING REMARKS 



The transition (0,3/2) <-> (1, 1/2) in the quantum-dot 
lithium is very similar to the Eermi-liquid - Wigner-solid 
transition in an infinite 2DESC3. This becomes evident if 
to plot results of quantum-dot and 2DES calculations in 
the form of an expansion in powers of Ti: energy in units 
e c i vs A -2 / 3 in dots, and energy in units e 1 JrT s vs r s 
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FIG. 1. Convergency of the energy of the lowest state with 
(Ltot, Stot) = (1,1/2). The energy unit is Tuvq. The curves 
are labeled by two numbers {E t h, N m b s ), where E t h is the 
threshold energy (in units huJo), and Nmbs is the number of 
many-body states involved in the expansion (0) . 
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FIG. 2. Low-lying energy levels with the total angular mo- 
mentum Ltot = 1 in a 3-electron quantum dot. The energy 
unit is hujQ. The five lowest levels are labeled by their total 
spin Stot and the degeneracy g. The insets schematically show 
the structure of the lowest-level wave function, at sm all (le ft 
inset) and large A (right inset), for details see Section HI C, 
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FIG. 4. Low-lying energy levels in a 3-electron quantum 
dot, for the total angular momentum L to t from to 2 and 
for all total spin states. Only the one lowest-energy level is 
shown for each pair of quantum numbers (Ltot, Stot)- The 
energy unit is Twjq. The levels are labeled by the pair of 
quantum numbers (Ltot, Stot) and the degeneracy g. 
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FIG. 3. Low-lying energy levels with the total angular mo- 
mentum Ltot = in a 3-electron quantum dot. The energy 
unit is huio- The five lowest levels are labeled by their total 
spin Stot and the degeneracy g. The inset schematically shows 
the stru cture of the lowest-level wave function, for details see 
Section |IIIC[ 




FIG. 5. 



Energy difference between the lowest states 
in a 3-electron quantum dot as a function 
of the interaction parameter A = Zo/as- The energy unit is 
Hljo- The transition occurs at A = A c ~ 4.343. 
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FIG. 6. Heat capacity of a 3-electron parabolic quantum 
dot, as a function of the temperature kT/huo and the inter- 
action parameter Zq/os. 




FIG. 7. Total electron density (spin-up plus spin-down) in 
the state (Ltot, Stot) = (0,3/2) at different A. Triangles show 
the positions of the classical radius (|l^) at corresponding val- 
ues of A. 






FIG. 8. Pair-correlation function Pjj of the state with 
Ltot = 0, Stot = 3/2, Sf* = +3/2, at A = 0.1, 2.0, 4.0, 
and 8.0, from up to down (|r'|/'o = Rd/h = 0.38, 1.04, 1.32, 
and 1.66, respectively). 
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FIG. 9. Total electron density (spin-up plus spin-down) in 
the state (1,1/2) at different A. Triangles show the positions 
of the classical radius ( |l9| ) at corresponding values of A. 




FIG. 10. Density of spin-up polarized electrons in the state 
(1,1/2) at different A. Triangles show the positions of the 
classical radius (h9|) at corresponding values of A. 
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FIG. 12. Pair-correlation function Pjj of the state with 
Liot = 1, Stot = 1/2, Sf* = +1/2, at A = 0.1, 2, 4, and 8, 
from up to down (|r'|/fo = Rd/h = 0.38, 1.04, 1.32, and 1.66, 
respectively) . 



FIG. 11. Density of spin-down polarized electrons in the 
state (1,1/2) at different A. Triangles show the positions of 
the classical radius ([L9|) at corresponding values of A. 
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FIG. 13. Pair-correlation lunction of the state with 
Ltot = 1, Stot = 1/2, S\ ot = +1/2, at A = 0.1, 2, 4, and 8, 
from up to down (|r'|// = R c i/h = 0.38, 1.04, 1.32, and 1.66, 
respectively) . 



FIG. 14. Pair-correlation function P-\y of the state with 
Ltot = 1, Stat = 1/2, S\ ot = +1/2, at A = 0.1, 2, 4, and 8, 
from up to down (|r'|// = R c i/h = 0.38, 1.04, 1.32, and 1.66, 
respectively) . 
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FIG. 15. Electron density in the states (1,1/2) (the ground 
state at A < A c ) and (0,3/2) (the ground state at A > A c ) at 
the transition point A = A c = 4.343. 
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FIG. 16. A "volume" -"pressure" diagram, A/Iq versus 
(10/A) 4 , for the ground state of a quantum dot lithium in the 
vicinity of the transition point A = A c = 4.343 (solid curve). 
Dashed curves show the "volume"- "pressure" diagrams for 
the states (0,3/2) (upper curve) and (1,1/2) (lower curve). 
Insets schematically show the distribution of electrons in the 
ground state on both sides of the transition point. In a real 
system the solid curve may have a hysteresis. 



TABLE I. Energies of the states (1,1/2) and (0,3/2) calcu- 
lated in this work (exact diagonalization) and in Ref.E3 (QMC, 
multilevel blocking algorithm). 



A Stot — 3/2 


Stot — 3/2 a 


Stot = 1/2 


Stot = l/2 a 


2 


8.3221 


8.37(1) 


8.1651 


8.16(3) 


4 


11.0527 


11.05(1) 


11.0422 


11.05(2) 


6 


13.4373 


13.43(1) 


13.4658 


no data 
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15.5938 


15.59(1) 


15.6334 


no data 


10 


17.5863 


17.60(1) 


17.6279 


no data 
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